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Multiple Scattered Radiation Emerging from Continental Haze Layers. 

1: Radiance, Polarization, and Neutral Points 

By George W. Kattawar, Gilbert N. Plass, and Stephen J. Hitzf elder 

ABSTRACT 

The complete radiation field including polarization is calculated by 
the matrix operator method for scattering layers of various optical thick- 
nesses. Results obtained for Rayleigh scattering are compared with those 
for scattering from a continental haze. Radiances calculated using Stokes 
vectors show differences as large as 23% compared to the approximate scalar 
theory of radiative transfer, while the same differences are only of the 
order of 0.1% for a continental haze phase function. The polarization 
of the reflected and transmitted radiation is given for a wide range of 
optical thicknesses of the scattering layer, for various solar zenith 
angles, and various surface albedos. Two entirely different types of 
neutral points occur for aerosol phase functions. Rayleigh-like neutral 
points (RNP) arise from the zero polarization in single scattering that 
occurs for all phase functions at scattering angles of 0° and 180°. For 
Rayleigh phase functions, the position of the RNP varies appreciably with 
the optical thickness of the scattering layer. At low solar elevations 
there may be four RNP. For a continental haze phase function the position 
of the RNP in the reflected radiation shows only a small variation with the 
optical thickness, and the RNP exists in the transmitted radiation only for 
extremely small optical thicknesses. Another type of neutral point (NRNP) 
exists for aerosol phase functions. It is associated with the zeros of the 
single scattered polarization which occur between the end points of the 
curve; these are called non-Rayleigh neutral points (NRNP) . There may 
be from zero to four of these neutral points associated with each zero 



2 


of the single scattering curve. They occur over a range of azimuthal 
angles, unlike the RNP which are in the principal plane only. The position 
of these neutral points is given as a function of solar angle and optical 
thickness. 


The authors are with the Physics Department, Texas A&M University, 
College Station, Texas 77843. 
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I. Introduction 

The radiation field produced by the photons which have undergone 
multiple scattering in a planetary atmosphere is in general polarized even 
though the incoming solar radiation is essentially unpolarized. Thus a 
complete description of the radiation scattered or transmitted by an atmos- 
pheric layer requires a specification of the degree and direction of polari- 
zation and the ellipticity in addition to the radiance. 

Elegant and elaborate mathematical solutions have been developed for 

1 2 3 4 

isotropic and Rayleigh scattering by Chandrasekhar, Sekera, ’ Kourganoff, 

and others. Numerical values for the various parameters of the diffuse 

radiation from a planetary atmosphere with Rayleigh scattering have been 

given in a set of tables published by Coulson et al.^ These tables are 

limited to optical thicknesses of unity or less. Numerical results for 

Rayleigh scattering have also been presented by Herman and Browning, Dave 

and Furukawa,^ Kahle,^’^ Howell and Jacobowitz,^*^ Fymat and Abhyankar,^^ 

12 

and Adams and Kattawar among others. Among these authors only Dave and 
7 9 

Furukawa and Kahle give radiance values for an optical thickness as large 

as 10. Most of the published results are for a surface albedo of zero. 

13 

More recently Plass ^ al have given radiance values out to large optical 

14 

thicknesses and for a range of values of the surface albedo. Plass et al 
have also calculated the radiance within optically deep absorbing media. 
Calculations of the polarization of the radiation are given only in Ref. 5, 

6, 7, 10, and 11 and only Dave and Furukawa^ give polarization values for 
an optical thickness as large as 10. 

The radiation derived by scattering for a haze layer has been calcu- 
lated by several different methods including Monte Carlo, iterative, doubling 
and matrix operator. Relatively few of these studies have included the 
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calculation of the polarization effects. Polarization values for the 

radiation scattered from various models of hazes have been reported by: 

t. ^ m 15,18,20,21,23 , j 16,17,19,22 

Kattawar and Plass and Plass and Kattawar 

24 25 26 

(Monte Carlo method) and by Dave, Herman et al, and Hansen (iter- 

24 

ative and doubling methods). Only Dave has given results for the 

ellipticity and direction of polarization in addition to the radiance 

^nd polarization values, but only for a size parameter of 10 and up to 

26 

an optical depth of 10. Hansen gives the polarization of the reflected 

light only over a complete range of optical depths . The Monte Carlo re- 
15-23 

suits cover a wide range of phase functions, optical depths, plane 

and spherical geometry, boundary surface including that of the ocean- 

atmosphere, but are averaged over a range of solid angle and necessarily 

have a certain statistical fluctuation in the results. 

Matrix operator theory is an entirely rigorous and practical method 

for the solution of the equations of radiative transfer. A brief history 

13 

of the method has been given by Plass ^ al . Results calculated by 

this method include the radiance of the reflected and transmitted light 

13 27 

for Rayleigh scattering layers and for maritime haze and the interior 

14 

radiances in optically deep absorbing media scattering by the Rayleigh 

28 

and a continental haze model. 

The matrix operator method has been generalized so that all four 
components of the Stokes vector are calculated. From this the radiance, 
degree and direction of polarization, and the ellipticity of the scattered 
radiation is obtained. The results are presented in the two parts of 
this paper for scattering from Rayleigh and continental haze layers. 

The radiance calculated from the approximate scalar theory is compared 
with the results obtained from the complete Stokes vector for both 
Rayleigh and haze scattering. The neutral points where the polarization 
is zero are calculated for both of these types of scattering. 
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II. Computational Aspects 

There is no difficulty in principle in generalizing the matrix opera- 

13 

tor method described by Plass ^ al to calculate the complete radiation 
field. The radiance is replaced by the four component Stokes vector in all 
equations, while the phase function is replaced by the fourth order phase 
matrix. All of the equations are formally the same with these and other 
appropriate substitutions. 

There are however, certain practical problems in the numerical solution 

of these equations. The storage requirements in the computer can become 

quite large and it is essential to minimize these as much as possible. 

These requirements can be greatly reduced by the use of the full symmetry 

relations which exist for the phase matrix of a spherical polydispersion as 

29 

presented by Kattawar e^ al. 

The order of the quadrature used is determined by the accuracy with 
which the phase function can be normalized. The order of the quadrature is 
considered acceptable if the normalization of the phase function differs 
from unity by less than one part in 10 . A symmetry preserving correction 
matrix is added so that normalization is achieved to the full accuracy of 
the computer. Great care was taken when adding the correction matrix. In 
extreme cases negative radiance values can be obtained even though the flux 
is conserved. 

A forty-second order Lobatto quadrature was used for the calculation 
presented in this article. Two different checks were made of these cal- 
culations against known results for scattering by the Rayleigh phase 
matrix. First the results for large optical depths were compared with the 
semi-infinite results given by Chandrasekhar K The Stokes vectors computed 
from the formulae on page 261 of Chandrasekhar agree with our matrix opera- 
tor results to one part in 10 , which is the accuracy of the functions 
tabulated by Chandrasekhar. Second, careful comparisons were made with 
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3 

the published tables of Coulson e^ ^ for a wide range of the parameters , 
Excellent agreement (to one or two places in the last tabulated figure) was 
obtained for every element of the Stokes vector except in a few regions 
where there are obviously small errors in the Coulson tables. For example, 
the parameter Q for the transmitted radiation is given in the tables as 
-0.00034 when x = 0.25, = 1, y“l, when the correct value is exactly 

zero . 

III. Radiance for Rayleigh Scattering Approximate Scaler vs. Stokes Results 

The upward radiance at the top of a plane parallel layer reflecting ac- 
cording to the Rayleigh phase matrix was calculated from the matrix operator 
theory and is shown in Fig. 1. The cosine of the solar zenith angle 
y^ = 1 and the surface albedo A = 0; the radiance is shown at various optical 
depths X from 0.0019 to 4,096. These results are presented here mainly for 
later comparisons with those calculated for haze layers and from the approximate 
scalar theory. This figure shows a change from limb brightening to limb dark- 
ening as the optical depth of the scattering layer increases. 

The diffuse downward radiance as observed at the bottom of the scattering 
layer is shown in Fig. 2 for y^ = 1. The maximum radiance is transmitted through 
such a scattering layer for optical thicknesses from about 1 to 4. The trans- 
mitted radiance becomes very small for the largest optical depth, t = 1,024, 
shown in the figure. 

The accurate radiance values calculated from the complete theory using the 
four-component Stokes vector are compared with the approximate values obtained 
from the scalar radiative transfer theory which uses only a single equation 
involving the radiance. The fractional difference in these two results 
AR/ R (R is exact radiance from Stokes vector; AR is. the exact radiance minus 
the approximate scalar radiance) is shown in Fig. 3 for the reflected radiation 
when y =1. The absolute value of the fractional difference between the two 
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results is small for small values of t. Increases to a maximum for t around 
unity, and then decreases to a limiting value as t increases further. The 
maximum absolute value of the fractional difference at x = 1 can be as large 
as 0.1 for a zenith angle of observation corresponding to V » 0.15. The 
maximum absolute value decreases to about 0.06 as t becomes large. 

The fractional difference in these two results is shown in Fig. 4 for the 
radiation transmitted through a layer of optical thickness t. The absolute 
value of 4R/R is a maximum for optical thicknesses around unity and then starts 
to decrease rapidly as x increases further and is quite small for x ^16. The 
largest difference that we found is when the sun is near the horizon and the 
observer is looking close to the horizon, when = 0.03785, P * 0.03785, 

X = 1, (|) = 0°, the exact radiance is 0.001625, while the scalar theory gives 
0.001998, for a difference of 23.0%. 

Radiance values calculated from the approximate scalar theory of radiative 
transfer for Rayleigh scattering can be in error by up to 23% at optical thick- 
nesses around unity and for the reflected radiation may be in error by 6 to 10% 
for all T ^1. Whether an error of this size is significant in a particular 
problem can only be determined by the proposed use of the results. 

IV. Radiance for Haze L Approximate Scaler vs. Stokes results 

Calculations were made of the complete radiation field scattered from the 

top of and transmitted through the bottom of a plane parallel continental haze 

31 

layer. The haze L model proposed by Deirmendjian is used. The number of 

2 1/2 

particles with a given radius is proportional to r exp (-15 . 1186r ), where r is 

the particle radius. The mode radius is 0.07 P®. A wavelength of 1 = 0.55 pm 

and real and imaginary parts of the index of refraction n^ = 1.55 and U 2 = 0.05 

were assumed as a reasonable approximation to the appropriate values for real 

continental hazes. The single scattering phase matrix was calculated by the 

32 29 

method described by Rattawar and Plass and Kattawar ad . The Mie scat- 
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tering matrix was obtained for size parameters x = 2 x r/A from 0.05 to 30 using 
a 75th order Gauss quadrature. The calculated single scattering albedo is 
0.7171. The single scattering function obtained from this phase matrix is shown 
in Fig. 5. There is a fairly strong maximum in the backward direction and a 
slight maximum at the rainbow angle. 

The upward radiance reflected from the top of a plane parallel haze L layer 
of optical thickness t is shown in Fig. 6 for * 0.85332 and A ■ 0. The 
azimuthal angle (j) is measured from the incident plane that contains the incident 
ray and the zenith direction. Here (j> = 0° on the left half of the figure and 
<j> = 180° on the right half; the antisolar point occurs in the right half of the 
figure. Thus as V varies from the left to the right in Fig. 6, we obtain the 
radiance values from the solar horizon to the nadir and then back through the 
antisolar point to the antisolar horizon. In all cases in this article the 
incoming solar flux is normalized to unity in a direction perpendicular to the 
solar beam. A maximum in the curve for the reflected radiance can be observed 
in Fig. 6 for all optical thicknesses near U = 0.853 and <f) = 180°. This is 
caused by the glory, the maximum in the single scattering curve for 180° scat- 
tering. The curve shown for T = 16 is a limiting curve which does not change 
further on the scale of the graph for larger optical thicknesses. 

The diffuse downward radiance transmitted through the bottom of a haze L 
layer of optical thickness t is shown in Fig. 7 for = 0.85332 (31.43°) 
and A = 0. As T varies from the left to the right in Fig. 7, we obtain the 
radiance values from the solar horizon through the solar point to the zenith and 
then back to the antisolar horizon. At very small optical thicknesses there is 
the expected maximum around the solar direction from the strong forward scat- 
tering of the phase function in addition to maxima at each horizon. As the 
optical thickness of the layer increases, the maximum around the solar direction 
noticably decreases and shifts toward the zenith when t is as large as 8, 

When T » 16, the radiance has a maximum at the zenith and no longer shows any 
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variation with azimuthal angle on the scale of this graph. 

In order to study the variation of the radiance with solar angle, the 
upward and downward radiance for = 0.3.3817 (79.15°) is shown in Figs. 8 and 
9. Tlie upward radiance for small optical thicknesses varies by four orders of 
magnitude between the solar horizon and the nadir and increases by two orders of 
magnitude from the nadir to the antisolar horizon. 

The downward radiance shown in Fig. 9 is particularly interesting as it 
Illustrates how the position of maximum radiance moves from the horizon to the 
zenith as the optical thickness of the haze layer increases. When t = 0.00195 
and 0.0625, the maximum is at the solar horizon. For larger optical thicknesses 
of the layer the angles at which the maximum radiance occurs is given in Table I 
for solar zenith angles of 31.43° and 79.15°. 

In the previous section it was shown that there may be differences as large 
as 23% in the radiance calculated from the approximate scalar theory compared to 
the correct value for the Rayleigh phase matrix. Are these differences as large 
for other types of phase functions? The haze L model was used for the calculation 
of the radiance by both the approximate scalar theory and by the correct Stokes 
vector theory. The results are compared in Table II. A comparison of a large 
number of these results for different values of t, and <t> showed extremely 

good agreement between the radiance results from these two theories. The 
largest differences found were of the order of 0.1%. The differences are much 
smaller than is the case for Rayleigh scattering, because of the smaller size of 
the off-diagonal matrix elements for the haze than for Rayleigh scattering. 

V. Polarization for Rayleigh Scattering 

In order to understand and appreciate the complexity of the polarization 
curves for haze L, it is necessary to study first the results for Rayleigh 
scattering. A few such results are given in this section. 
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The polarization of the photons reflected from a Rayleigh layer is shown in 
Fig. 10 for = 0.18817 (79.15°), A = 0, and <f> = 0°. The polarization of 
single scattered photons from a Rayleigh scattering function is well known to be 
zero for scattering angles of 0° and 180° and to be a maximum at 90°. The curve 
for T = 0 in Fig. 10 refers to single scattered photons and is zero when 
h = 0.18817 (®-= 79.15°, the antisolar direction; scattering angle = 180° and 
is unity when V = 0.98214 (9 = 10.85°; scattering angle = 90°). As soon as 
there is a finite optical thickness, the single neutral point at the antisolar 
direction splits into two neutral points, one on each side. As the optical 
thickness, increases the neutral points move farther away from the original 
position. One neutral point eventually reaches the antisolar horizon and then 
reappears just above the solar horizon. The polarization of the reflected 
radiation in general decreases as increases in regions away from the neutral 
points (a broad region around the nadir in the example shown in Fig. 10), but 
may increase appreciably as "T Increases in a region around the antisolar direction. 
The curve shown for x = 4 is a limiting curve that does not change on the scale 
of this figure as x increases further. 

The neutral points that develop for Rayleigh scattering have often been 
discussed in the literature. The most complete curves previously published 
showing the position of the neutral points as a function of the optical thick- 
ness of the layer are those of Dave and Furukawa^. However, their curves are 
for the transmitted radiation only and for a single solar zenith angle (84.26°). 
They pointed out for the first time that there are four neutral points over a 
certain range of optical thickness near unity at this particular solar zenith 
angle . 

The neutral points for the reflected radiation as a function of the optical 

thickness of the reflecting layer are shown in Fig. 11 forh ^ = 0.85332 (31.43°), 

y = 0.11333 (83.49°), andy = 0.03785 (87.83°). The ordinate of these figures 
o o 

is the cosine of the nadir angle of observation and the abscissa is the optical 
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thickness of the reflecting layer. The Babinet point (between the antisolar 
direction and the nadir) moves toward the nadir as the optical thickness in- 
creases. It reaches a limiting position at an optical thickness of about 10 that 
does not change further as the optical thickness further increases . Curves are 
shown for several surface albedos (A = 0 and 1) . The position of the Babinet 
point is relatively insensitive to A. 

The Brewster point (between the antisolar direction and the antisolar 
horizon) moves from the antisolar point toward the antisolar horizon as the 
optical thickness increases. When the svin is sufficiently near the horizon, the 
Brewster point touches t^ie antisolar horizon at some value of the optical thick- 
ness. As the optical thickness further increases, the neutral point reappears 
on the solar horizon and is now called the Arago point (between the solar hori- 
zon and the nadir) . A further increase in the optical thickness causes the 
Arago point to move up from the solar horizon. The position of the Brewster and 
Arago points depends slightly on the surface albedo A. When the position of the 

Brewster and Arago points is plotted as in Fig. 11, the curve is continuous 

across the horizon. 

The position of the neutral points when ■ 0.53786 (57.46°) and “ 0.18817 
(79.15°) is given in Fig. 12. The position of the Brewster point for 6^ * 57.46° 
is relatively sensitive to the surface albedo A. 

The polarization of the transmitted photons for = 0.18817 (79.15°), A = 0, 
and (j) - 0° and 180° is given in Fig. 13. In this case the single scattered 
photons have zero polari ation when y = 0.18817 (9 = 79.15°, the solar direction; 
Scattering angle = 0°). As the optical thickness of the layer increases, the 
two neutral points move out from y = 0.18817. Eventually one neutral point 
appears on the antisolar horizon. When 1 < x <2.5, there are four neutral 
points, and, when x > 2.5, there are none. The polarization curve for x = 4 has 

no zeros and has relative maxima near the zenith and at both horizons. This 

curve is not quite S3rmmetric around the zenith. 
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The shape of the polarization curve for radiation transmitted through a 
layer of optical thickness t Is entirely different when t = 16. This is an 
asymptotic curve which is the same for all layers with larger optical thicknesses. 
The polarization is zero at the zenith. After the photon has passed through the 
diffusion region in the scattering layer, it has lost all memory of its original 
direction; thus the polarization at the zenith of the radiation leaving the 
lower surface of the medium must be zero, since there can be no preferred direc- 
tion. Furthermore the curve must be independent of the azimuthal angle. It is 
perhaps somewhat surprising that the radiation emerging near the horizon from a 
very thick Rayleigh scattering layer has a polarization of the order of 10%. 

The position of the neutral points of the radiation transmitted through a 
Rayleigh layer of given optical thickness is shown in Figs, l^f - 17 for four 

different solar angles. The position of the neutral points ji^hen = 0.85332 

/ 

(31.43°) is given in Fig. 14, a case typical of those with/the sun near the 

j 

zenith. The Babinet point (between the solar direction ayhd the zenith) moves 

/ 

from the solar direction toward the zenith as the optical thickness increases, 

I 

reaches a maximum displacement at about t = 1.3, and ^then moves back toward the 

/ 

solar direction to meet the Brewster point at t = 5/. 2. There are no neutral 
points at larger values of t when A = 0. The Brewster point (between the solar 
direction and the solar horizon) moves from the solar direction toward the solar 
horizon as t increases, reaches a maximum displacement at about = 2.1, and 
then moves back toward the solar direction, to meet the Babinet point. There is 
a considerable difference in the curves for a surface albedo A = 0 and A = 1, 
but little difference between A = 0 and A = 0.2. There are never more than two 
neutral points for this solar angle. 

The same curves, but for ] 1 ^ = 0.53786 (57.46^) are given in Fig. 15. The 
unusual aspect of this example is that only Babinet and Brewster points, exist 
when the surface albedo A = 0, while an Arago point (between antisolar horizon 
and zenith) exists when 0.68 <t> 3.4 and A = 1. 
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The neutral points when = 0.18817 (79.15^) are shown in Fig. 16. When 
the sun is this close to the horizon an interesting new phenomenon appears (as 
was pointed out for the case = 0.1 by Dave and Furukawa^). Four neutral 
points (two Arago, one Brewster, and one Babinet) appear over the range of 
optical thickness from 1.2 to 1.8 when A = 0. There is another range of optical 
thickness from about 2^0 to 2.5 where there are two Babinet points and no others. 
The curves when = 0.03785 (87.83^} are given in Fig. 17. There is one 
Babinet and one Brewster point for A = 0 when t< 0.0015; one Babinet and one 
Arago point when 0.0015 < t < 0.96; one Babinet, one Brewster, and two Arago 
points when 0.96 < r < 1.1; two Babinet and two Arago points when 1.1 < x < 1.6; 
two Babinet points when 1.6< x < 1.7; no neutral points when x > 1.7. 


VI. Polarization for Haze L Scattering 

The polarization of the radiation scattered by aerosol layers exhibits 
various new features which are not found in scattering from Rayleigh type parti- 
cles. One major difference is that the polarization of single scattered radi- 
ation from Rayleigh particles is zero only at the end points (scattering angle 
(j) = 0^ and 180^). In general, when there is moderate or no absorption, the 
single scattered polarization of radiation scattered from aerosols exhibits one 
or more zeros in between the end points of the curve. 

The polarization of single scattered photons from haze L is shown in 
Fig. 18 as a function of the cosine of the scattering angle 0 , Four different 
curves show the polarization on succesively magnified scales. The following 



points should be noted: 1. the polarization reaches a maximum value of about 

^ 0.53 at 0 '*158^; 2. the polarization is zero for scattering angles of 3.53^ 

H ■, : .. / 

i andf25.6^; 3. the polarization is very small (less than 10 for 0*^ < 0 <3.53°; 

the absolute value of the polarization is small (less than 0.0064) for 3.53 

■ ' ■■ \' . 

< 0 < 25 . 6 °. 
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The single scattered photons have zero polarization whenever the SGattering 
angle in three dimensions is either 3.53° or 25.6°. There are two curves for 
each solar zenith angle in Fig. 19, which give the corresponding polar coordi- 
nates of these particular scattering angles. Any point on a curve gives a value 
for V (cosine of the zenith angle) and of ^ (azimuthal angle) for which the 
polarization is zero. There are two curves for each solar zenith angle; the 
small curve near the left hand edge of the figure is for the scattering angle 
3.53°, while the other curve is for 25.6°. When the sun is near the zenith the 
polarization can be zero for © = 3.53° and 25.6° only for the radiation trans- 
mitted through the bottom of the layer, but not for the reflected radiation; on 
the other hand, when the sun is close to the horizon, there are in addition 
angles at which the polarization is zero for the radiation reflected from the 
top of the layer. 

These neutral points will be called NRNP («on-Rayleigh neutral points) , 
since they have no counterpart in Rayleigh scattering. The NRNP arise from 
zeros in the single scattering polarization curve at scattering angles between 
0° and 180°. On the other hand the RNP (Rayleigh neutral points) arise from the 
zeros that always occur in the single scattering polarization curve at the end 
points (0=0° and 180°). As our examples will show, the behavior of the polari- 
zation near NRNP and RNP is quite different as the optical thickness increases . 

The polarization of the radiation reflected from and transmitted through 
haze L layers is Shown in the following examples. The polarization of the 
reflected photons is given in Fig. 20 when = 0.85332 (31.43°). The curves in 
the upper part of the figure are for the principal plane (^ =0° and 180°), 
while those in the lower part are for (|» = 30° and 150°. For single scattering 
(t = 0) there is only a RNP in the principal plane when - 0.85332. There is 
no NRNP for the reflected radiation at such a high sun angle. As t increases 
two neutral points develop around the single scattering RNP, just as occurs for 
scattering from Rayleigh type particles. In the present case the two neutral 
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points never move as far from = 0.85332 as they do in the case of Rayleigh 

scattering. They move less than 3° as is shown in Table III. 

The largest values of the polarization in Fig. 20 occur in two regions on 
either side of the RNP near angles of observation that correspond to a single 
scattering angle of 158^ from the incident solar beam. 

The position of the RNP as a function of the optical thickness of the 
reflecting haze L layer is shown in Fig. 21 for five different solar zenith 
angles. In all cases there is a relatively small change (usually less than 5^) 
in the position of the RNP as the optical thickness increases from zero. The 
small variation in these neutral points is in sharp contrast to the much larger 
changes that occur for scattering from a Rayleigh phase function. 

The polarization of the reflected photons for y^ = 0.85332 in the plane at 
right angles to the principal plane ((fr = 90°) is given in the top part of Fig. 22 
(with the polarization of the photons transmitted through the layer in the lower 
part of the figure) . 

An example of the polarization of the reflected photons from a haze L layer 
when the sun is near the horizon is given in Fig. 23 for the case y^ == 0.188166 
(79.15°). The RNP exist near the antisolar direction (y = 0.188166, (|) = 180°). 

As the optical thickness increases these move away from the antisolar direction 
as shown in Fig. 21 and in Table III. In addition, at this low solar angle, 
there is also a single NRNP near y = 0.26 (75.2°; ^ = 0°). As previously dis- 
cussed and shown in Fig. 19, the scattering angle for the single scattered 
reflected photons is 25.6° in this direction; this scattering angle corresponds 
to a zero in the single scattering function for haze L and causes this NRNP. 

The position of this NRNP changes only slightly as the optical thickness of the 
reflecting layer increases; Table III shows that it changes from 75.2° for x = 0 
to 75.9° for a large optical thickness. 
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The polarization of the photons transmitted through a haze L layer is shown 
in Fig. 24 when = 0.85332 (31.43°) and ^ = 0° and 180°. First there is a RNP 
at V = 0.85332, <ji = 0° for single scattered photons (t = 0). This RNP no longer 
exists for radiation transmitted through layers of even very small optical 
thickness, presumably because of the extremely small values for the single 
scattered polarization at small scattering angles. 

In addition there are two NRNP which are near p = 0.544 (57.0°) and 0.9952 
(5.6°). As the optical thickness increases, these NRNP change their position as 
indicated in Fig, 24 and Table III. They finally disappear for an optical 
thickness between 8 and 16. 

The polarization curves for the plane ^ - 30° and 150° are shown in Fig. 25^ 
For small optical thicknesses there are two NRNP near P = 0.66 and 0.99 at the 
angles indicated in Fig. 19. The position of one of these NRNP changes appreciably 
when T = 4 and both NRNP disappear for an optical thickness between 4 and 8. 

The NRNP occur over a range of azimuthal angles, while the RNP never occur 
outside of the principal plane. 

The nature of the polarization curves is quite different at low solar 
elevations as shown in Fig. 26 for the case P^ = 0.188166. There is a RNP at 
P = 0.188166 and cf) = 0° for single scattered photons; it disappears for layers 
with very small optical thicknesses. There is one NRNP at P = 0.593 (53.6°) 
for single scattered photons, the position indicated in Fig. 19. As the optical 
thickness of the haze layer increases, this NRNP has moved only slightly to 
57.6° for T= 2 (see Fig. 26 and Table III). Interestingly a second neutral 
point appears at the horizon when t = 1,4 and there are two NRNP when 1.4 < x< 2.3. 
When T = 2 the two NRNP are at 57.6° and 70.1° (see Fig. 26 and Table III). 

There are no neutral points for the transmitted radiation when'^ > 2.3. 

The position of the NRNP for five solar angles is given in Fig, 27. The 
curves vary greatly in shape depending on the solar angle. There are one or two 
NRNP depending on the optical thickness. There is a suggestion in the curve for 
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= 0.53786 that four NRNP exist over a very narrow range of optical thickness. 
VII. Conclusions 

The polarization of the radiation reflected from and transmitted through a 
scattering layer is compared for Rayleigh and haze L phase matrices. In both 
cases there are neutral points arising from the zeros of the polarization of 
single scattered photons at scattering angles of 0° and 180°. The angular 
position of these Rayleigh-like neutral points (RNP) in the sky shows appreci- 
able variation with the optical thickness of the scattering layer for a Rayleigh 

phase matrix, but only a small variation for haze L. In fact RNP exist for the 
transmitted radiation through a haze L layer only when the optical thickness is 
extremely small. RNP exist only in the principal plane. 

Another type of neutral point exists for non-Rayleigh phase functions which 
is associated with the zeros of the polarization for single scattering which 

occur between the end points of the curve; these are called non-Rayleigh neutral 

points. For a given optical thickness there may be zero, one, two, or even four 
of these neutral points associated with each of the zeros of the polarization 
for single scattering. NRNP occur over a range of azimuthal angles starting 
from the principal plane. 

A comparison of radiances calculated from the complete theory of radiative 
transfer using Stokes vectors with those obtained from the scalar theory shows 
that differences of the order of 23% may be obtained for Rayleigh scattering, 
while the largest difference found for a haze L phase function was of the order 
of 0.1%. The scalar theory can apparently be used with confidence to calculate 
radiances for hazes and aerosols. 
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Table I. 

Angles for maximum radiance and polarization. 


Solar 

Zenith 

Angle 

Optical 

thickness 

Zenith, angle for 
maximum downward 
radiance 
(<}>=0°) 

Maximum Polari- 
zation of down- 
ward radiance 
((j)=180o) 

Maximum Polari- 
zation of upward 
radiance 
((j)=180O) 




Zenith Angle 

Nadir 

angle 

Nadir 

angle 

31.43° 

0 



53.56° 

10.36° 


0.00195 

32.22° 


53.56° 

10.36° 


0.00781 

32.22° 


53.56° 

10.36° 


0.0625 

32.22° 


53.62° 

10.23° 


0.25 

32.11° 


53.77° 

9.80° 


1 

32.00° 


54.48° 

8.54° 


2 

31.68° 


55.18° 

7.73° 


4 

30.80° 


56.01° 

7.65° 


8 

25.58° 


56.15° 

7.65° 


16 

4.36° 


56.15° 

7.61° 

79.15° 

0 

90.00° 

78.87° 

57.25° 



0.00195 

90.00° 

78.81° 

57.25° 



0.00781 

90.00° 

78,64° 

57.18° 



0.0625 

90.00° 

77.52° 

56.63° 



0.25 

80.87° 

75.26° 

55.60° 



1 

73.50° 

68.78° 

54.27° 



2 

60.59° 

55.60° 

54.41° 



4 

37.81° 

46.68° 

54.48° 



8 

11.87° 

46.45° 

54.48° 



16 

2.69° 

90.00° 

54.48° 
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Table II. 

Comparison, of exact and approximate scalar- theory radiance values. 


T 


P 


Downwelling 
or upwelling 
radiation 

Scalar 

Theory 

Stokes Vector 
Theory 

Percent 

Difference 

0.25 

1 

1 

- 

Up 

2.8621 -3 

2.8631 -3 

0.03% 

0.25 

1 

0.03785 

- 

Up 

1.5092 -2 

1.5089 -2 

-0.02% 

0.25 

0.03785 

0 . 89034 

60° 

Ip 

9.6381 -4 

9.6367 -4 

-0.01% 

0.25 

0.11333 

0.99574 

90° 

Down 

2.0010 -3 

2.0006 -3 

-0.02% 

1, 

0.53786 

0.89034 

60° 

Up 

7.9171 -3 

7.9146 -3 

-0.03% 

1 

0.85332 

0.11333 

150° 

Down 

1.3534 -2 

1.3528 -2 

-0.04% 

2 

0.18816 

1 

30° 

Up 

3.7611 -3 

3.7596 -3 

-0.04% 

2 

0.11333 

0.03785 

90° 

Down 

6.3774 -4 

6.3706 -4 

-0 . 11% 

4 

0.53786 

0.98574 

60° 

Up 

9.7842 -3 

9.7811 -3 

-0.03% 

4 

0.85332 

0.40452 

180° 

Down 

9.6320 -3 

9.6283 -3 

-0.04% 

8 

0 . 18816 

1 

0° 

Up 

3.9351 -3 

3.9336 -3 

-O.v'4% 

8 

0.03785 

0.11333 

30° 

Down 

6.8123 -6 

6.8172 -6 

0.07% 

16 

0.53786 

0.99574 

120° 

Up 

9.2782 -3 

9.2769 -3 

-0.01% 

16 

0.18816 

0.03785 

150° 

Down 

5.9750 -7 

5.9793 -7 

0.07% 
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Table III. 


_ Optical 
thickness 

Solar 

zenith 

angle 

Nadir angle of 
neutral points of 
reflected photons 

Zenith angle of 
neutral points of 
transmitted photons 

m » 


4>=180° 

<|)=180° 

4)=0° 

(j)=0° 

jb. 00195 

31.43° 

31.4° 

31.4° 

57.0° 

5.6° 

0.00781 


31.4° 

31.4° 

57.0° 

5.4° 

1 p.0625 


31.6° 

31.3° 

57.0° 

5.3° 

M it 

0.25 


31.8° 

31.2° 

57.1° 

5.7° 



32.5° 

30.5° 

57.5° 

6.5° 

^'2 


33.6° 

29.3° 

57.8° 

7.2° 

T’* 


34.0° 

28.8° 

57.9° 

9.6° 

ii 


34.0° 

28.8° 

43.4° 

28.1° 

f6 


34.0° 

28.8° 

- 

- 



(}i = 0° 

(j)=180° 

(|)=0° 

(j)=0° 

1 J. 00195 

79.15° 

75.2° 

78.8°; 79.6° 

53.6° 

-*• 

’ 6.00781 


75.2° 

75.5°; 80.7° 

53,6° 

- 

iib.0625 


75.0° 

76.7°; 82.0° 

53.8° 

- 

i ■ 

*0.25 


75.2° 

75.6°; 83.4° 

54.2° 

- 

p “« 


75.8° 

75.1°; 83.9° 

54.6° 

- 

1 : 

2 } 


75.9° 

75.2°; 83.7° 

57.6° 

70.1^ 

^ 4 


75.9° 

75.2°; 83.7° 

- 

- 

LI 


75.9° 

75.2°; 83,7° 

- 

- 

16 


75.9° 

75.2°; 83.7° 

- 

- 


/I ^ 

t : 


I ^ 
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Fig. 1. 


Fig. 2. 


Fig. 3. 


Fig. 4. 
Fig. 5. 
Fig. 6. 


Fig. 7. 


Captions for Figures 

Upward radiance reflected from a plane parallel layer with Rayleigh 
phase matrix as a function of the cosine of the nadir angle y. The 
cosine of the solar zenith angle y^ = 1; the surface albedo A = 0. 
Curves are given for layers with various values of the optica] 
thickness x. The incoming solar flux is normalized to unity in a 
direction perpendicular to the solar beam. 

Downward radiance transmitted through a plane parallel layer with 

Rayleigh scattering as a function of the cosine of the zenith angle 

y. y = 1; A = 0. 
o 

AR/R (where R is exact radiance from Stokes vector; AR is the exact 
radiance minus the approximate scalar radiance) for the reflected 
light from a Rayleigh scattering layer for y^ = 1, A = 0. 

AR/R (see caption to Fig. 3) for the transmitted light through a 
Rayleigh scattering layer for y^ = 1, A = 0. 

Single scattering function obtained from phase matrix for haze L 
(continental haze) as a function of the scattering angle. 

Upward radiance reflected from a plane parallel haze L layer as a 
function of the cosine of the nadir angle y for - 0.85332, 

A = 0. Tl\e left hand portion of the graph is for an azimuthal 
angle (}> = 0° aAd the right hand portion for <j> = 180° , The solar 
horizon is at the left of the figure, the nadir at the center fol- 
lowed by the antisolar point and the antisolar horizon at the right. 
Curves are given for layers with various values of the optical 
thickness x. 

Downward radiance transmitted through a haze L layer as a function 
of the cosine of the zenith angle y for y^ = 0.85332, A * 0. The 
azimuthal angle ifi = 0° for the left hand portion of the graph and 
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Fig. 

Fig. 

Fig. 

Fig. 


Fig. 

Fig. 

Fig. 


180° for the right. The solar horizon is at the left of the figure, 
followed by the solar point, the zenith, and the antisolar horizon 
on the right. 

8. Upward radiance for haze L, p = 0.18817, A = 0, and <)) = 0° and 

o 

180°. (See caption to Fig. 6). 

9. Downward radiance for haze L, = 0.18817, A = 0, and (|) = 0° and 
180°. (See caption to Fig. 7). 

10. Polarization of the upward radiation for Rayleigh scattering for 

U = 0.18817, A = 0, and (j) = 0° and 180° for layers of various 
o 

optical thicknesses. (See caption to Fig. 6). 

11. Position of Babinet, Brewstsr, and Arago neutral points for upward 

radiation reflected from a Rayleigh scattering layer when 

V =0.85332, 0.111333, and 0.03785 and A = 0 and 1. The ordinate 
o 

is the cosine of the nadir angle and the abscissa is the optical 
thickness of the reflecting layer. The Brewster point appears be- 
tween the antisolar direction and the antisolar horizon and the 
Arago point appears between the solar horizon and the nadir. 

12. Position of Babinet, Brewster, and Arago neutral points for upward 

radiation reflected from a Rayleigh scattering layer when 

y = 0.53786 and 0.18817 and A = 0, 0.2, and 1. (See caption to 
o 

Fig. 11). 

13. Polarization of the downward radiation transmitted through Rayleigh 
scattering layers of various optical thicknesses for y^ = 0.18817, 

A = 0, and (|) = 0° and 180°. (See caption to Fig. 7). 

14. Position of Babinet and Brewster neutral points for downward radi- 
ation transmitted through a Rayleigh scattering layer when 

y = 0.85332 and A = 0, 0.2, and 1. The ordinate is the cosine of 
0 

the zenith angle and the abscissa is the optical thickness of the 


transmitting layer. 
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Fig. 15. 


Fig. 16. 


Fig. 17. 


Fig. 18. 


Fig. 19. 


Fig. 20. 


Fig. 21. 


Position of Babinet, Brewster, and Arago neutral points for down- 
ward radiation transmitted through a Rayleigh scattering layer 

when u = 0.53786 and A = 0, 0.2, and 1. The Brewster point appears 
o 

between the solar direction and the solar horizon while the Arago 
point appears between the antisolar horizon and the zenith. The 
ordinate is the cosine of the zenith angle and the abscissa is the 
optical thickness of the transmitting layer. 

Position of Babinet, Brewster, and Arago neutral points for down- 
ward radiation transmitted through a Rayleigh scattering layer when 

V = 0.18817 and A = 0 and 1. (See caption to Fig. 15). 
o 

Position of Babinet, Brewster, and Arago neutral points for down- 
ward radiation transmitted through a Rayleigh scattering layer when 

y = 0.03785 and A = 0 and 1. (See caption to Fig. 15). 
o 

Polarization of single scattered photons from haze L. The two in- 
sets show the polarization for scattering angles near 0°. The 
abscissa is the cosine of the scattering angle. 

Angles at which the polarization of single scattered photons from 
haze L is zero. Each curve is for a particular value of the cosine 
of the solar zenith angle y^. The ordinate is the cosine of the 
zenith angle y and the abscissa is the azimuthal angle (|). The 
upper part of the figure is for photons reflected into an upward 
direction, while the lower part is for photons transmitted through 
the scattering layer in a downward direction. 

Polarization of the upward radiation for haze L, y^ = 0.85332, 

A = 0, and (ji = 0° and 180° (upper curves) and (ji = 30° and 150° 
(lower curves) for layers of various optical thicknesses. (See 
caption to Fig. 6). 

Position of Rayleigh-like neutral points for photons reflected from 
a haze L layer. Curves are given for five different solar zenith 



- 27 - 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


angles. The ordinate is the cosine of the nadir angle and the 
abiscissa is the optical thickness of the reflecting layer. 

22. Polarization of the upward and downward radiation for haze L for 

= 0.85332, A = 0, and (j)'= 90° for layers of various optical 
thicknesses. (See captions to Figs. 6 and 7). 

23. Polarization of the upward radiation for haze L, = 0.188166, 

A = 0, and (j> = 0° and 180° for layers of various optical thicknesses. 
(See caption to Fig. 6). 

24. Polarization of the downward radiation for haze L, u = 0.85332, 

A = 0, and (f) = 0° and 180° for layers of various optical thicknesses. 
(See caption to Fig. 7). 

25. Polarization of the downward radiation for haze L for = 0.85332, 

A = 0, and ^ = 30° and 150° for layers of various optical thicknesses. 
(See caption to Fig. 7). 

26. Polarization of the downward radiation for haze L for = 0.188166, 

A = 0, and = 0° and 180° for layers of various optical thicknesses. 
(See caption to Fig. 7). 

27. Position of non-Rayleigh-like neutral points for photons transmitted 
through a haze L layer. Curves are given for six different solar 
zenith angles. The ordinate is the cosine of the zenith angle and 
the abscissa is the optical thickness of the scattering layer. 
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